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Abstract 

The MiniBooNE experiment at Fermilab has amassed the largest sample to date of 
7r°s produced in neutral current (NC) neutrino- nucleus interactions at low energy. 
This paper reports a measurement of the momentum distribution of 7r°s produced 
in mineral oil (CH2) and the first observation of coherent ir° production below 
2 GeV. In the forward direction, the yield of events observed above the expectation 
for resonant production is attributed primarily to coherent production off carbon, 
but may also include a small contribution from diffractive production on hydrogen. 
Integrated over the MiniBooNE neutrino flux, the sum of the NC coherent and 
diffractive modes is found to be (19.5 ±1.1 (stat) ±2.5 (sys))% of all exclusive 
NC 7T° production at MiniBooNE. These measurements are of immediate utility 
because they quantify an important background to MiniBooNE's search for v„ — * v e 
oscillations. 
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1 Introduction 

Neutral current (NC) n° production is the single largest ^-induced back- 
ground to neutrino experiments measuring —>■ v e oscillations in the E v ~ 
1 GeV range, including the search recently performed by the MiniBooNE ex- 
periment pQ . NC 7T° events can mimic v e signal events when, for example, one 
of the two photons associated with the 7r° — > 77 decay is not detected. This 
can happen when a photon exits the detector before showering or does not 
have enough energy to initiate a shower. Estimating the rate of such back- 
grounds relies on knowledge of neutrino induced NC tt° production at low 
energy (E u < 2 GeV). 
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Pion production from the scattering of low energy neutrinos on nuclei prin- 
cipally occurs through two mechanisms. The larger contribution comes from 
incoherent processes where the neutrino interacts with one of the nucleons 
in the nucleus. In the MiniBooNE energy range this mainly consists of the 
excitation and subsequent pionic decay of baryonic resonances (such as the 
A(1232)). Additionally there is a small but non-negligible yield from coherent 
scattering where the neutrino interacts with the entire nucleus leaving it in its 
ground state. Because of the necessarily small momentum transfer, coherent 7r° 
events are more forward peaked than their resonantly produced counterparts. 

To predict the full spectrum of it production reliably, it is important to char- 
acterize the resonant ir° contribution, as it is the dominant source of ir°s. Mini- 
BooNE models resonantly produced NC 7r° events using the Rein and Sehgal 
(RS) model [2J as implemented in version 3 of the Nuance event generator [3] 
assuming an iV — > A dipole form factor with axial mass, M r £ s = 1.1 GeV/c 2 . 
For MiniBooNE, 95% of resonant NC tt° production is predicted to occur 
via the A(1232), but seventeen higher mass resonances with their interfer- 
ences also contribute in the model. 85% of the resonant NC 7r° production at 
MiniBooNE should occur on carbon with the remaining 15% on hydrogen. To 
predict both the kinematics and yield of coherently produced 7r° events, Mini- 
BooNE uses the RS coherent production model [I], implemented in Nuance 
with the relevant axial mass set to Mjf = 1.03 GeV/c 2 . The model predicts 
coherent 7r° production to be 30% of the total NC exclusive n° production 
in MiniBooNE. The Nuance implementation differs from the RS model in 
two important ways. First, resonances are decayed isotropically, which is not 
strictly correct. Events are reweighted to match the RS model based on the 
A decay angle in its rest frame with respect to its momentum vector. Second, 
Rein and Sehgal describe an absorption factor, which scales the coherent pro- 
duction cross section for NC ir°s, while in Nuance, absorption is implemented 
as part of the final state interaction (FSI) model. The designations of resonant 
and coherent are set prior to any FSI, which means that rescattered events 
with a 7T° in the final state may be misclassified in Nuance, as would be the 
case when a coherently produced 7r° rescatters elastically through a resonance. 

Calculating 7r° production cross section in either case becomes complicated for 
several reasons. In the case of resonant tt° production, the neutrino-nucleon 
cross section requires knowledge of the appropriate transition form factors. Us- 
ing the CVC hypothesis [5|6] , vector form factors can be reliably inferred from 
electron scattering data; however, axial- vector form factors are not well known 
and rely heavily on the use of PCAC [7]. Added to this, for neutrino-nucleus 
scattering, one must additionally calculate the propagation of 7r°'s through the 
target nucleus as this can change both the identity of the pion as well as its 
kinematics. At resonance energies, the pion-nucleon cross section is large and 
hence the pion produced in the resonance decay has a non-negligible proba- 
bility to re-interact before exiting the target nucleus. This must be properly 
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accounted for in any useful simulation, because the experimental observable 
is the 7T° only after it has exited the target nucleus. 

In the case of coherent scattering, the situation is even more ambiguous. Calcu- 
lations of coherent scattering cross sections have been performed using detailed 
neutrino-nucleon resonance production models and subsequent hadron-hadron 
interactions [8119] . Such a recent calculation [10] finds a value of 14% for the 
ratio of coherent to incoherent scattering for the NC process investigated here. 
Alternatively, one can circumvent some of the complexities of these dynamics 
by invoking Adler's PCAC theorem [TT], which directly identifies the coher- 
ent scattering cross section with the elastic scattering of pions on the same 
nucleus. This procedure [4"f8|12lll3|ll4|ll5j works well at higher energies, but ap- 
pears to fail at low energies where PCAC-based calculations typically predict 
a substantially larger fraction of coherent scattering. For example, using data 
on 7r + C scattering [16] . one would infer a coherent fraction of roughly 50%. 
In addition to the differing theoretical approaches and large range in coherent 
scattering predictions, the K2K experiment has recently, and somewhat sur- 
prisingly, reported no evidence for charged current coherent 7r + production at 
1.3 GeV [17] . For these reasons, experimental measurements of coherent pion 
production are critical to our understanding of this complex process, and are 
especially important at low energy. 

To date, there are only a few measurements of NC tt° production in the 1- 
2 GeV energy range, conducted on a variety of different targets materials, 
which together consist of about 3000 events [18 ; 19 ; 20 : 21]. NC coherent 7r° 
production data are even more sparse. Although cross section measurements 
exist at higher energies [22], there are no measurements of NC coherent tt° 
production below E u < 2 GeV, which is an important region for neutrino 
oscillation experiments. 

The remainder of this paper describes the MiniBooNE experimental setup, the 
identification and reconstruction of NC ir° events, a measurement of the overall 
yield of NC 7r° production in mineral oil as a function of 7r° momentum, and a 
direct measurement of the coherent 7r° production fraction in this data sample. 
Understanding coherent production is critical in reproducing the observed 
angular spectrum of these events. Together, this technique and resultant 7r° 
constraints provide important input to the MiniBooNE oscillation analysis [1] 
and update previously reported work on this subject [23, 241125] . 



2 The Experiment 

The MiniBooNE z/ M beam results from the decays of secondary particles (mostly 
pions) that are produced by interactions of 8 GeV protons from Fermilab's 



4 



Booster incident on a beryllium target. The detector, 541 m downstream of 
the beryllium target, is a 12.2 m diameter spherical tank filled with 800 tons 
of pure mineral oil. The tank is separated into two regions: an inner volume 
with a radius of 575 cm, and a 35 cm thick outer volume. An optical bar- 
rier provides the separation of the two regions, and also serves as the support 
structure for 1280 equally-spaced 8-inch photomultiplier tubes (PMTs) that 
give 10% photocathode coverage of the inner volume. PMT hits have a thresh- 
old of ~0.1 photoelectrons, and are recorded in a 19.2 fis window around every 
1.6 fis neutrino beam spill. An additional 240 8-inch PMTs mounted in the 
outer volume act as a veto shield to detect charged particles entering or exiting 
the detector. Three meters of dirt above the detector give a 60% reduction 
in cosmic ray flux, and with appropriate selection cuts, the outer veto region 
rejects more than 99.95% of the cosmic rays observed in the detector. 

The energy spectrum peaks at 700 MeV and extends to approximately 
3 GeV. Integrated over the neutrino flux, approximately 7% of the neutrino 
interactions in MiniBooNE are predicted to be NC exclusive 7r° production. 



3 Event Selection and Analysis 

The selection of NC ir° events begins with a set of simple pre-reconstruction 
cuts (or pre-cuts) that exactly match those used in the MiniBooNE v e ap- 
pearance analysis pQ. Each candidate must have only a primary event, with 
no evidence of a secondary event consistent with a muon decay electron. This 
eliminates the vast majority of CC interactions. Each event is also required 
to have more than 200 PMT hits in the main tank, well above the maximum 
number of hits observed at the muon-decay endpoint. Each event must have 
fewer than 6 hits in the veto region, which eliminates cosmic rays and neutrino 
events not contained in, or originating outside of, the detector. Additionally, 
all events must be in the 1.6 fis beam spill window, although, after the pre- 
ceding cuts are applied, almost no events exist outside the beam window. 

Each event that passes the pre-cuts is then reconstructed under three hypothe- 
ses |2l)] : muon, electron, and it . The reconstruction is based on the expected 
distribution of Cerenkov and scintillation light in the mineral oil. Under the 
muon and electron hypotheses, the times and charges of hit PMTs are used to 
reconstruct a single track fitting the location and time of the neutrino interac- 
tion as well as the energy and direction of the charged lepton track. Electron 
and muon tracks are distinguished by the energy deposition per unit length 
and by the sharpness of the Cerenkov ring edge (electron rings are fuzzier than 
muon rings due to scattering and the formation of electromagnetic showers). 
The 7T° hypothesis requires a two track fit which fits the location and time 
of the neutrino interaction as well as the energies, directions, and conversion 
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Cut Level 



Survival 



Pre-cuts 



100% 



r <500 cm cut 



79.4% 



log(£ e /£ M ) > 0.05 

log(Ar/£e) > 

80< m 77 <200 MeV/c 2 



39.5% 



51.5% 



50.1% 



Table 1 



NC 7T° event survival fraction (in percent) for each selection requirement relative 
to the pre-cuts. The cuts are applied progressively from the top of the table down. 

distances of the two photons. Each photon track fit assumes that the light is 
distributed in the manner of an electron track. Each fitted hypothesis produces 
a likelihood (£ M , £ e , £ w ), and the logs of the ratio of likelihoods from different 
hypotheses are used for particle identification. The 7r° fit is run in two ways: 
with a floating invariant mass, m 77 , and with m 77 fixed to the 7r° mass. The 
7T° parameters, with the exception of the invariant mass, are obtained from 
the fixed mass fit as this provides the most accurate estimate of the true ir° 
kinematics. 

The reconstructed parameters and likelihoods allow for further selection. Be- 
cause their decay photons shower like electrons, 7r° events should look more like 
electrons than muons, and overall these events should look more like a 7r° than 
an electron. Therefore, particle identification cuts requiring log(£ e /£ M ) > 0.05 
and log(£ w /£ e ) > are applied. Additionally, a fiducial volume cut requires 
the reconstructed position of the event to be within 500 cm of the detector 
center. These cuts produce a very clean sample of 7r° events with a signal to 
background ratio of ~30. With an additional selection on the invariant mass 
(80 < m 77 < 200 MeV/c 2 ), the ir° efficiency predicted by the Monte Carlo 
(MC) is 39.5%, as detailed in Table E After all cuts, the MiniBooNE data 
set consists of 28,000 NC 7r° events produced in 5.6xl0 20 protons on target, 
which is the largest sample of NC 7r° events yet collected at these energies. 

Once selected, the 7r° candidate events are divided into bins of reconstructed 
7T° momentum and the MC is used to unsmear the data (i.e., to reverse the 
effects of momentum resolution and inefficiency and thus obtain the "true" ir° 
production rates as a function of momentum). A matrix is formed by divid- 
ing MC events into bins of true momentum versus reconstructed momentum 
and counting true 7r° events in each bin. Events in this matrix must pass all 
selection cuts including the mass window cut. A MC event is defined to be 
a 7t° event if at least one decayed ir° exists in the final state. This definition 
includes both NC and CC events, although most CC events are eliminated by 
the pre-cuts which exclude events with electrons from muon-decay. The event 
count in each bin is divided by the total number of tc° events in that true 
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Fig. 1. Top: Results of the ir° unsmearing in bins of momentum. The dark points 
show the unsmeared data 7r° momentum distribution and the light points show the 
uncorrected MC 7r° momentum distribution. The unsmeared data error bars contain 
all sources of error propagated through the unsmearing, while the MC error bars 
result solely from finite MC statistics. Bottom: The reweighting function, formed 
by taking the ratio of the two distributions in the top plot (data/MC). 



momentum bin (including events that did not pass the reconstruction cuts, 
which are included in the denominator to correct for the cut efficiency). This 
matrix is inverted to form the unsmearing matrix. Next, a vector is formed in 
the data by separating the 7r° candidate events into the same reconstructed 
bins. The background rate in each data bin is estimated using the MC and 
the data event yields are scaled to remove this estimated background. The 
product of the unsmearing matrix and the data vector is the unsmeared data 
vector. While in many applications this kind of matrix unsmearing can be 
unstable, leading to large uncertainties in the unsmeared quantities [27], in 
this case, with the initial matrix largely diagonal, the process is quite reliable. 
This assertion is supported both by MC closure tests and by the reasonable 
size of the errors on the unsmeared data as propagated through the matrix. 

Figure [T] shows a comparison of the initial MC prediction to this unsmeared 
data distribution. The ratio of the two distributions forms a reweighting func- 
tion which is used to scale MC 7r° events as a function of true momentum. 
This distribution reflects the extent to which the starting MC does or does 
not predict the measured momentum spectrum of 7r° events in MiniBooNE. 
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By construction, this momentum reweighting fixes the discrepancy between 
data and MC in reconstructed 7r° momentum. Additionally, it improves agree- 
ment in many key kinematic distributions. Figure [5] shows relatively normal- 
ized data to MC comparisons for both the initial and corrected MC. The 
kinematic distributions shown are the cosine of the 77 opening angle, the 
photon energies, and the 7r° momentum. All distributions show marked im- 
provement (this is almost a tautology for the momentum distribution, except 
that it is binned more finely than the correction function and is in terms of 
reconstructed momentum). This reweighting of the MC is neither profound 
nor forbidden; it merely addresses a range of imperfections in the simulation, 
which may include contributions from the neutrino flux to the ir° production 
model. 



4 Results 



The 7T° candidate events in the momentum reweighted MC are divided into 
three templates: one each for resonant (res), coherent (coh), and background 
(bg) events. The resonant and coherent templates contain all true exclusive 
7T° event 9^1. In the case of the coherent template, this includes diffractive 
scattering off hydrogen [28J, which for MiniBooNE is predicted to be 16% of 
the coherent template. The background template consists of all other events, 
including some events that contain one or more decayed ir° not produced in 
the resonant or coherent channels. 

The templates are formed as a two dimensional (2D) distribution of E^l — 
cos8 n ) versus invariant mass m 77 . Use of this 2D distribution helps to break 
the degeneracy between the coherent and resonant templates in m 77 , and 
between the coherent and background templates in E n (l — cos9 n ). The angle 
9-n is defined to be the lab angle of the reconstructed ir° momentum vector with 
respect to the neutrino beam directional The 2D template binning is defined 
by dividing the ID distributions (-£^(1 — cos9 n ) and m 77 ) from the MC into 
variable-width bins of approximately equal numbers of events. Each fit has 
three parameters (x res , x coh and x bg ) which scale the template distributions 



By our definition, a true exclusive ir event must be generated by Nuance in 
either the coherent or resonant modes and there must be a decayed 7r° in the final 
state. 

2 The simpler angular function, cos 9 n , was tried in place of E n (1 — cos 9 n ) and found 
to have slightly poorer performance. This is attributed to the fact that E n (l — cos n ) 
has a more consistent shape across all ir° momenta for coherent events. 
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Fig. 2. Relatively normalized comparison of the uncorrected (dotted) and momen- 
tum corrected (solid gray) Monte Carlo to data (points with statistical error bars) 
for various 7r° kinematic distributions: a) the opening angle between the two pho- 
tons in the ir° decay, b) energy of the more energetic photon, c) energy of the less 
energetic photon, and d) the tt° momentum. The marked improvement in these 
kinematic distributions shows that the initial data-to-MC differences can largely be 
attributed to the original disagreement in 7r° production as a function of momentum. 

independently. The fit minimizes the following \ 2 '- 

r fdata / fres r^res \ fcoh^cohi fbg bg\]2 

2 _ Ui KJi X ~T~Ji X ^Ji X )\ Q\ 

i {°t ata ) 2 + (°r) 2 + (vt oh ) 2 + (<h 9 ) 2 



where ff is the fraction of total events of type a (where a is data, res, coh, 
and bg) in the ith bin and a? is the statistical uncertainty on that fraction. 
Since the fit is to the fractional distributions, it is a shape only fit and the 
sum of the fit parameters should be very close to unity, as is the case for all 
fits. 
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Fig. 3. Comparison of data and Monte Carlo in a) m 77 invariant mass, and b) 
E^il — cos 6 W ) after the coherent fraction fit. The resonant, coherent and background 
components are shown scaled by their fit parameters. The full MC fit is the sum of 
the three components. 

The fit is repeated for several different binnings. The number of bins in each 
ID projection is varied independently from 15 to 25, for a total of 121 different 
binning combinations. The final fit parameters are formed from the average 
of the parameters from the 121 fits. To determine the best overall production 
parametrization, the momentum correction and coherent fit are iterated, using 
the results of one as a correction to the inputs of the other. This procedure 
converges after only two iterations. Figure [3] shows the final fit plotted in the 
1 — cos^) projections. The fit coherent fraction is defined as: 



m 77 and E. 



coh 



X 



X' 

coh 



coh 



+ X r 



x 100%. 



(2) 



The fit finds F co h to be (19.5 ± l.l(stat)%. The MiniBooNE data clearly favor 
the presence of a coherent scattering component. The average confidence level 
(C.L.) of the fit is 7.14%, while the C.L. obtained when the coherent fraction 
is fixed to zero (x coh = 0) is 10~ 18 . The effects of the momentum reweighting 
are small, but not insignificant. If the momentum reweighting is not done the 
fit coherent fraction is 18.5%. 

One should note that the reported coherent fraction is specific to the Mini- 
BooNE neutrino spectrum and includes scattering off both carbon and hydro- 
gen nuclei in the mineral oil target. It has also been measured in the context 
of the RS-based Nuance generator [4][5] . with the aforementioned modifica- 
tion to the A decay angular distribution. This widely-used model predicts a 
coherent fraction of 30% for MiniBooNE. Of course, more recent calculations 
of coherent production [9|12|ll4j do predict a range of lower coherent fraction 
values for MiniBooNE. Figure H] compares the measured coherent fraction to 
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Fig. 4. The coherent fraction in the Rein and Seghal based MC vs. neutrino energy 
compared to this measurement. The solid line includes only carbon interactions, 
while the dotted line includes scattering off hydrogen with diffractive events counted 
as part of the coherent. The measured value is shown with error bars which indicate 
the total error on the measurement (vertical) and the spread in the participating 
neutrino energy distribution (horizontal). The shaded distribution is MC energy 
for neutrinos which produce NC 7r° events in MiniBooNE with arbitrary normal- 
ization. The coherent fraction predicted by Nuance integrated over all energies in 
MiniBooNE is 30%. 

the RS/Nuance prediction as a function of neutrino energy. The plot shows 
two predictions: one with both carbon and hydrogen scattering (dashed) and 
another which includes only carbon interactions (solid). The effect of hydrogen 
scattering is small compared to the precision of the measured coherent frac- 
tion. Using the MC to correct to a pure carbon target would yield a measured 
coherent fraction of (20.3 ± 2.8(stat))%. The shaded distribution shows the 
predicted neutrino energy spectrum for neutrinos which participate in NC n° 
production in MiniBooNE. 



5 Systematic Uncertainties 

Systematic uncertainties on the coherent fraction include choice of binning, 
background composition, momentum reweighting, neutrino flux, choice of anal- 
ysis cuts, and detector modeling. The binning systematic is deduced from the 
RMS on the coherent fraction from the fits to the 121 different binnings. To 
determine the background shape uncertainty, background events are divided 
into several classes and the production cross section of each class is randomly 
varied according to a Gaussian distribution (with standard deviations from 1 
to 40% depending on the estimate of the uncertainty in each process class). 
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This is repeated 5000 times and the background shape from each combination 
of variations is used in the template fit. The background uncertainty is given by 
the RMS of these 5000 fits. The reweighting error is determined by randomly 
varying the momentum reweighting function according to its errors, paying 
careful attention to bin-to-bin correlations. Estimating this uncertainty pro- 
ceeds as in the background case with 5000 random combinations and fits. The 
flux uncertainty results from varying parameters in the beam simulation. The 
analysis cuts error is determined by varying the cut point on the reconstructed 
variables. Finally, the detector model uncertainty is determined by fitting, as 
fake data, 70 data-sized MC samples which were simulated with random, but 
properly correlated variations in several detector response parameters. Since 
the 70 samples are statistically independent, the detector model error is: 

&det. model = \] RM S 2 - {(J fit ) 2 (3) 



where RMS is the root mean square of the 70 fits and ((Tfit) is the average fit 
error of 1.1%. 

Table [2] lists the uncertainties estimated from each source. The dominant 
source of systematic error is the detector model, which is largely due to the un- 
certainty in the reconstructed energy which is strongly correlated with several 
of the varied parameters. 



6 Conclusions 



Using a high statistics sample of events, MiniBooNE has measured the rate 
of NC 7T° production in mineral oil as a function of momentum and extracted 
a correction to the predicted production rate for this process. MiniBooNE 



Error Source Coherent Fraction (%) 

Binning 0.21 

Background Shape 0.64 

Reweighting 0.51 

Flux 0.06 

Analysis Cuts 0.51 

Detector Model 2.34 

Total Systematic Error 2.54 

Table 2 

Contributions to the systematic uncertainty in the MiniBooNE measured coherent 
fraction. 
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also reports the observed rate of coherent ir° production relative to the total 
exclusive 7r° production in the context of the RS model [I] as implemented in 
Nuance[3]. The coherent fraction is found to be (19.5±l.l(stat)±2.5(sys))% 
for the MiniBooNE flux and target. This should be compared to the 30% 
fraction predicted by the RS-basedNuANCE model, a value significantly higher 
than the measurement reported here. The fit to the MiniBooNE NC 7r° sample 
excluded, with high confidence, the possibility of no coherent contribution to 
7T° production at MiniBooNE energies. 

In the MiniBooNE u e appearance oscillation analysis [I], both the ir° momen- 
tum correction and the measured coherent fraction were used to reweight the 
MC for a more accurate estimation of tc° misidentification in MiniBooNE. By 
correcting the ir° production with data as described in this paper, the error 
on the overall 7r° production (an input to the oscillation analysis) goes from 
~25% (the quadratic sum of a ~20% flux error and a ~20% NC tt° cross 
section error) down to 5%. This represents an important improvement in the 
sensitivity to v e appearance. 
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